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Abstract 

We consider a nonlocal reaction-diffusion equation as a model for a population structured 
by a space variable and a phenotypical trait. To sustain the possibility of invasion in the case 
where an underlying principal eigenvalue is negative, we investigate the existence of travelling 
wave solutions. We identify a minimal speed c* > 0, and prove the existence of waves when 
c > c* and the non existence when < c < c*. 
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1 Introduction 

1.1 Setting of the problem 

In this paper we are interested in propagation phenomena for nonlocal reaction-diffusion equations 
of the form 



d t n(t,x,y)-A Xty n(t,x,y) = lr(y - Bx ■ e) - J k(y - Bx ■ e,y' - Bx ■ e)n(t,x,y') dy' J n(t,x,y), 

(1) 

where (x,y) € R d x K, e e S"* -1 , B > 0, r : K -> R and k : M 2 -> R+. 

Such equations have appeared in some population dynamic models, see [38], [41], [40], [35]. 
In this context, n(t, x, y) denotes a density of population structured by a spatial variable x € K d 
and by a phenotypical trait y6M. The population is then submitted to four essential processes: 
spatial dispersion, mutations, growth and competition. The spatial dispersion and the mutations 
are modelled by diffusion operators. The growth rate of the population at a location x and trait 
y is given — for all times t — by r(y — Bx ■ e), where r is typically negative outside a bounded 
interval. This corresponds to a population living in an environmental cline: to survive at the 
location x, an individual must have a trait close to y = Bx ■ e. Therefore, to be able to invade 
the environment, the population needs to evolve. Finally, we consider a logistic regulation of the 
population density that is local in the spatial variable and nonlocal in the trait. In other words, we 
consider that there exists an intra-specific competition (for e.g. food) at each location, which may 



1 I3M, Universite de Montpellier 2, CC051, Place Eugene Bataillon, 34095 Montpellier Cedex 5, France. E-mail: 
malfaro@math.univ-montp2.fr 

2 Equipe BIOSP, INRA Avignon, Domaine Saint Paul, Site Agroparc, 84914 Avignon Cedex 9, France. E-mail: 
jerome.coville@avignon.inra.fr 

3 Centre d'Ecologie Fonctionnelle et Evolutive, UMR 5175, CNRS, 1919 Route de Mende, 34293 Montpellier, 
France. E-mail: raoul@cefe.cnrs.fr . 



1 



depend on the traits of the competitors. For a rigorous derivation of this model from individual 
based models, we refer to [16]. In Section 2, we will discuss in more details the biological aspects 
of our work. 

The existence of global solutions for the Cauchy problem and of non trivial steady states for 
(1) have been investigated respectively in [41] and [2]. Also, numerical simulations (see e.g. [38], 
[41], [40]) show that the population can either go extinct, or propagate while adapting to local 
environments. 

The aim of this work is to analyze this propagation phenomena through the study of travelling 
front solutions. The travelling front solutions are particular solution of (1) describing the transition 
at a constant speed c from one stationary solution to another one. Such solutions have proved in 
numerous situations their utility in describing the dynamics of a population modelled by a reaction 
diffusion equation. In the case of the classical Fisher-KPP equation 

d t n - An = (1 - n) n, (2) 

we refer among others to [23], [33], [3] [43]: there exists planar fronts <fi(x.e — ct) connecting to 
1, for all speed c > c* = Moreover, the minimal speed of the front c* corresponds to the so 
called spreading speed of propagation. Travelling front solutions in heterogeneous versions of (2) 
with periodicity in space, in time, or more general media are studied in [12], [28], [44], [6], [8], 
[36], [37]. Nonlocal versions of (2) where the Laplace operator is replaced by a nonlocal operator 
are studied in [19], [17, 18], [42]. For very general reaction diffusion equations, we refer to [7] for 
a definition of generalized transition waves and their properties. 

It is worth noticing that when the competition term is replaced by a local (in x and y) density 
regulation, equation (1) becomes the following heterogeneous reaction diffusion equation 

d t n{t, x, y) - A x>y n(t, x, y) = (r(y - Bx ■ e) - h(y - Bx ■ e)n(t, x, y)) n(t, x, y), (3) 

which was recently investigated by Berestycki and Chapuisat [5]: they prove the existence of a 
critical speed c* , for which there exists a travelling front of (3) for any speed c > c*. 

As far as nonlocal equations of the form (1) are concerned, far less seems to be known in 
the literature. The travelling wave analysis has been done either using a formal Hamilton-Jacobi 
approach [15] for a model close to (1), or for a population structured by one variable only, that is 
n = n(t,x), submitted to a nonlocal competition [9]. To our knowledge, there is no result on the 
existence of travelling waves for (1) and related models. 



1.2 Assumptions and main results 

As suggested by the numerical simulations mentioned above, we expect that during an invasion, 
the population adapts locally to the environmental gradient. To observe travelling waves of (1), 
we therefore perform the change of variable 

h(t, x, z) = n(t, x, z + Bx ■ e). (4) 

Then (1) is recast as 

dtfi(t, x, z) — £ (n)(t, x, z) = I r(z) — / k(z,z')h(t, x, z')dz' ) h(t, x, z), (5) 



where £(n) := A x h+(B 2 + 1) n zz -2Bd z (V x n-e). Since (1) is invariant under any rotation in R , 
without loss of generality we can assume that e = e\ so that £ (n) = A x h+ (£> 2 + l) n zz — 2Bd XlZ n. 
This operator is elliptic, since the associated matrix has only positive eigenvalues. Looking after 
travelling wave solutions, we search a speed c and a profile u(x, z) such that h(t, x, z) := u(x ■ 
e\ — ct, z) = u(xi — ct, z) solves (5). For convenience, we drop the numerical subscript and write 
x instead of x\. Hence, we are looking after (c, u{x, z)) such that 

— £(u)(x, z) — cu x (x, z) — I r(z) — / k(z, z')u(x, z 1 ) dz' I u(x, z) in M 2 , (6) 
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where 



£(u) := u xx + (B 2 + l)u zz - 2Bu x 
Throughout the paper, we make the following assumption. 



0,0/ 



Assumption 1.1 (Structure of r and k) Function r is in the Holder space C loc 
< 9 < 1, and there is 5 > such that 

1 



for some 
(7) 



Vz e R, r(z) < ^ - 6z 2 . 

o 

Function k is in the Holder space C"f c (R 2 ) and th ere are k > 0, k + > such that 

V(z,z')eR 2 , k~ <k(z,z')<k + . 

Let us next introduce a principal eigenvalue problem that is necessary to enunciate our main 
result. For more details on principal eigenvalue problems in general domains we refer to [13], [14] 
and the references therein. 

Definition 1.2 (Principal eigenvalue problem) We denote by (A^r^) € R x C°°(R) the 
solution of the principal eigenvalue problem 



- (B 2 + 1) A,r°, (z) - r(z)r°, (z) = A ^ (z) for all z e 
r^(z)>0 forallzeR, T%(0) = 1. 



(8) 



Observe that in the case where r(z) = 1 — Az 2 , A > 0, we have A^, = y A (_B 2 + 1) — 1 and 



r^ c (z) = exp ^— y -q^pj'yJ is a Gaussian profile. 

We first state that as soon as A^ > 0, extinction of the population occurs. 
Proposition 1.3 (Extinction) Assume A^ > 0. For any initial population n° such that 

n°(x,y) 



< oo, 



T^iy-Bx) 

the solution of (1) with initial condition n° goes extinct exponentially fast as t — > oo. ■ 

n(t,x,y) 



L°°(R 2 ) 



0(e 



-\'t\ 



Next, we state our main result: as soon as A^, < 0, invasion waves exist. Precisely, the following 
holds. 



Theorem 1.4 (Travelling waves) Assume A^ < and define 

c* := 2 



B 2 + l' 

Then the following holds. 

(i) For all c> c* , there exists a positive u € C 2 (R 2 ) solution of 



(9) 



with 



£{u)(x, z) — cu x (x, z) = I r(z) — / k(z, z')u(x, z') dz' I u(x, z) in R 2 , (10) 



(11) 



I/ l{(x,z)e(-oo,o]x[-i/,i/]} < u(x,z) < Ce 



-Kz 
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for some v > 0, C > 0, K > 0, and 



\\u(x, OIU ->x->-+oo 0, / u(x, z) dz ^x^ +00 0. 
Additionally, when c > c* , i/iere exists /x < smc/i £/ia£ 



(12) 



u( V-B 2 + l xH — J z ] „ 



(13) 



(m) When < c < c*, i/iere is no positive solution of (10) smc/i £/ia£ lim inf x ^+oo u(x, 0) = and 
m(x, z) < V'( z ) / or some ?/> G L 1 (R). 

1.3 Comments 

On the extinction case. The proof of Proposition 1.3 is elementary and we now give the proof. 
The result is a consequence of the parabolic comparison principle satisfied by the local equation 



Indeed, one can check that Me A °°*r^ (?/ — Bx.e) and n(t, x, y) are respectively a super- and a 
sub-solution of (14) with ordered initial data (for M large enough). 

On the construction of waves. Let us first comment on a major difficulty in the construction of 
travelling fronts. When the competition term is replaced by a local (in x and y) density regulation, 
many techniques based on the comparison principle — such as some monotone iterative schemes 
or the sliding method [11] — can be used to get a priori bounds, existence and monotonicity 
properties of the solution. Since integro-differential equations with a nonlocal competition term 
do not satisfy the comparison principle, it is unlikely that such techniques apply here. 

It turns out that the considered problem here has some similarities with the case of a population 
structured by a spatial variable only, that is n = n(t,x), submitted to a nonlocal competition as 
studied in [9] (see also [1]). In this work, the construction of a travelling front is based on a 
sequence of approximating problems on intervals (— a n ,a n ), with a n — > oo. Due to the lack of 
comparison principle for the approximated problem, the construction of a solution is based on a 
topological degree argument, a method introduced initially in [10]. 

To construct our fronts, we adopt a similar strategy and consider a sequence of problems in 
growing boxes (— a„,a n ) x (—b n ,b n ), with a normalization at the origin. In order to make this 
strategy possible, a key point is to establish a priori estimates, independent on the size of the 
boxes, on the profile u, the speed c and, in particular, the tails of u when z is large. Due to the 
nature of the considered kernels here, an uniform estimate on u is obtained using a local pointwise 
L p estimate, whereas the uniform control on c is obtained by showing that our problem does not 
have a solution if the speed c is too large or if c = 0. Notice that the latter analysis may also be 
used to simplify the proof in [9] . 

Let us highlight that, in contrast with [10] and [9], it is far from obvious that the constructed 
travelling fronts are monotone w.r.t. x for x > large enough. Therefore, we shall need an extra 
work to catch the behavior (12) ass^ +oo. 

Notice also that the comprehension of the behavior of the wave as x —> — oo is quite involved. 
In the related case of the nonlocal Fisher-KPP equation, the positive steady state u = 1 may 
present, for some kernels, a Turing instability (see e.g. [24], [9], [1]). Such a situation may also 
occur in our context. 

Let us also mention that, although we construct fronts without relying on any monotonic 
properties of the profiles, it is suspected, as in the case of the nonlocal Fisher-KPP equation [9], 
that there exist monotone and non monotone travelling fronts. The understanding of such issues 
is quite challenging. 

Organization of the paper. In Section 2 we briefly describe the biological context of (1) and 
give an interpretation of our results. Then, we prove Theorem 1.4 in Sections 3-5. In Section 3, 



d t <t>{t,x,y) - A Xi y(j)(t,x,y) 



r(y - Bx.e)(j){t,x,y). 



(14) 
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we start by deriving some a priori bounds and then, using a Leray-Schauder topological degree 
argument, we construct a solution in a bounded box. In Section 4, we let the box tend to R 2 and 
obtain a wave, which turns out to be the one with minimal speed c — c* . We also show the non 
existence of waves with speed < c < c*. Lastly, we construct faster waves c > c* in Section 5. 



2 Biological interpretation of the results 

In this section we briefly precise the biological context of (1). 

In the present paper, we are interested in biological invasions involving darwinian evolution. 
Species invading new territories often face environmental gradients of e.g. temperature, luminosity, 
antibiotic chemicals. Experimentally, it is well documented that invasive species then evolve during 
their range expansion [22], [30], to adapt to local conditions. To understand the speed, or even 
the success of an invasion, one should thus consider the dispersion, birth and death processes, 
but should also take into account evolution [26], [31], [27]. Those questions become especially 
important in the context of the global warming [20], [21]: the favorable environmental conditions 
of many species move towards the north, implying important changes in species' range. It is also 
of great importance for the evolution of resistance of bacteria to antibiotics [27] . 

More generally, many evolutionary biology questions involve spatially structured populations, 
while most existing models either neglect the spacial structure of the population, or largely sim- 
plify it. New theoretical tools are then needed, and structured population models are natural 
candidates: they enable the modelling of all the biological phenomena mentioned above, and nu- 
merical simulations show that they are able to reproduce interesting features. Analyzing this type 
of model is however challenging, even in a homogeneous setting, see e.g. [29], [34]. This work, 
as well as the results of [15] are first steps in the mathematical understanding of the dynamics of 
those models. 

The main application of our result concerns asexual populations living in an environmental 
cline. The simplest model then writes 

0~ 2 0~ 2 

d t n(t, x, y) — £A x n(t, x, y) - -y A y n(t, x, y) = 

"max - -Jn-iV- hx f - -F? I x : V') d V) x > V), 

where o~ x , a m describe respectively the diffusion rate and the mutation rate of the population, 
is the strength of the selection, b is the steepness of the environment cline, and K the carrying 
capacity of the environment. After the rescaling n(t,x,y) = n(r max t, ^ 2 ^"' ax x, ^ 2 J'"° X y), we see 
that h solves (1) with 

9 



r(y) = l-Ay\ A=^y, B = f-b, 



The population then gets extinct if A (B 2 + l) > 1, while if A (B 2 + l) < 1, invasion fronts exist, 
with a minimal propagation speed (in the original variables) 



v 7 ^ 



max u x 




-1/2 



Remark 2.1 There exists thus only two dynamics: either the population gets extinct, or it succeeds 
to invade the whole territory. The situation of asexual populations is then very different from the 
case of sexual populations (see [32], [35]), where populations surviving with a limited range only 
are possible. 

Notice also that during invasions, the dispersion of individuals can evolve (see e.g. [39]). Our 
result does not apply to this problem, and we refer to [4-], [15] for such a situation. 
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3 The problem in a bounded box 
3.1 On some principal eigenvalue problems 

We first introduce some principal eigenvalue problems, whose eigenfunctions will serve as boundary 
conditions when stating the travelling wave problem in a bounded box. 

For v G [0, 8), where 6 > is as in Assumption 1.1, we denote by (A^,, rj^) the solution of the 
principal eigenvalue problem 

i-(B 2 + l)A z r^(z)~(r(z) + ^z 2 )r^(z)^X^(z) for all z G R 
|r^(z)>0 forallzGM, 1^(0) = 1. 

Notice that this definition is coherent with (8), and for any v G [0,6), we have A^, < A^, < 0. 
Also, for v G [0, 8) and 6 > 0, we define (A^ , T%) as the solution of the principal eigenvalue problem 

'— (B 2 + 1) A z T» b (z) - (r(z) + vz 2 ) r%(z) = \ v b T%(z) for all z G (-b,b) 

< m±b) = o (i6) 

T%(z) > for all z G (-6, b), T^(0) = 1. 

Let us observe that b < b' implies X^ < X b , < X b , and that X b — > X'^ as b — > oo. To construct 
the travelling waves, we will use the eigenfunctions r b , for b > as a boundary value. To 
bound from above those functions independently of b > 0, we will also use the functions T^ 3 . 
Notice that — max ze i r(z) < X^ 3 < xH 3 < 0. To show that T^ 3 is integrable, we define 

C := 
and 



C := max[_j !? ] T^ 3 , where z := and tp(z) := Cexp (— \J B i +1 z 2 Jk ) ' so ^( = '- z ) = ^ 



- [B 2 + 1) A^(z) - (r(z) + yz 2 ) Hz) > (Jz 2 - V(^) > 0, 

for all z G (— oo, —z) U (z, oo). Since in (-co, —z) U (z, oo) we have r(z) + ^z 2 < j — |z 2 < 0, the 

comparison principle then applies to (15) on (— oo, —z) U (z, oo), and yields T^ 3 (z) < ip(z) for all 
z G (— oo, —z) U (z, oo). As a result, for some constant which we denote again by C , we have 

^(z) < Cexp ^-*—J^j forallzGM, (17) 

which implies in turn that T^ 3 G L 1 (R). 

For a given 6o > 0, we now use a similar argument to control the functions T b (z) uniformly 

f /3(a£ /3 -a^ 3 ) I » / s/a A 1 
w.r.t. b G [6o, oo]. When z lies in [— z, z], where z := max < y - — 2— 3 1 \ 4s \ \ + ^ ) (1 the 

coefficients of the equations in (15) and (16) are uniformly bounded w.r.t. b G [60,00]. Therefore 
the Harnack inequality implies that there is C > such that T b (z) < C, for all z G [— z, z], all 

b G [60,00]. By the definition of z we see that, on the one hand, T 2 ^ is a super-solution for (15) 
and (16) — with v = | — on (— 00, — z) U (z, 00) and that, on the other hand, the comparison 
principle applies. Therefore, there exists C > such that 

T S b /3 (z)<CT 2 J 1 / 3 (z)<C\\r 2 ^ 3 \\ 00 for all z G M, feG [60,00]. (18) 

In particular, we have 

I T S b ' 3 (z)dz < m T := C [ T^ /3 (z) dz < 00 for all 6 G [6 , 00] . (19) 
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3.2 The problem in a box 

For a > 0, b > and e G (0, 1), we consider the problem of finding a speed c <G K and a real 
function u{x,z), defined for (x, z) G [—a, a] x [—6,6], such that 



— £(u)(x, z) — cu x (x, z) 



P{a,b,e) < 



l {u(x,z)>0 



} r(z) 



k(z 7 z')u(x, z') dz' u(x,z) 



u(x,z) = l {x= ^ a} (x)T S b /3 (z) 
u(0,0) = e, 

where Q := (—a, a) x (—6, b). The elliptic operator is given by 



in Q 
on dQ 



-£(u) = -u xx - (B 2 



2Bu, 



(20) 



If (c, u) is a solution achieving a negative minimum at (x m ,z m ) then, from the boundary 
conditions we deduce that (x m , z m ) lies in the interior of the rectangle, and that — £ (u) — cu x = 
on a neighborhood of (x m , z m ). The maximum principle thus implies u = u(x m , z m ), which cannot 
be. Therefore any solution of P(a, 6, e) satisfies u > and, by the strong maximum principle, 

u>0 and — £ (u)(x, z) — cu x (x, z) = (r(z) — J k(z, z')u(x, z 1 ) dz'\ u(x, z) in Q . (21) 



In the following of the section, we shall construct a solution to P{a, b, e) via a Leray-Schauder 
topological degree argument. To make this possible, we consider a family of problems as follows. 
For a > 0, b > and r G [0, 1] , we consider the problem of finding a speed c G K and a nonnegative 
real function u(x, z) such that 



' — £(u)(x, z) — cu x (x, z) 



P T (a,b) { 



= \ r(z) - t I k(z, z')u(x, z') dz' — 7(1 — t)m(.t, z) I u(x, z) in 



, u(x,z) = l {x= _ a} (x)T S b /3 (z) 



on dQ, 



where 7 > will be specified later (see Lemma 3.7). Note that Po(a, b) reduces to a local prob- 
lem, and that solving P(a, b, e) is equivalent to solving Pi (a, b) with the additional normalization 
condition u(0, 0) = s. 

Remark 3.1 A first natural idea to define a family of problems would be to consider —£{u) — cu x = 

t (r{z) — J^ h fc(z, z')u(x, z') dz'^j u. But then we cannot get a uniform w.r.t. < r < 1 control of 

the tails of u (see Lemma 3.4), which is crucial to derive e.g. a lower bound on the standing waves 
(see Lemma 3.6). This is the reason why we consider the family P T (a,b) as above. Therefore, the 
topological degree argument (see subsection 3.6) is rather involved and requires to analyze a whole 
family of local problems (see Lemma 3.8). 



3.3 A priori estimates for u 

We provide a priori bounds for the profile u of solutions to P T (a, b). When < r < 1/2, the local 
part of the equation shall be enough to derive Lemma 3.3. On the other hand, when 1/2 < r < 1, 
the nonlocal part is quite relevant and we first need to control the vertical mass of u, namely 

,6 

m(x) := / u(x,z)dz. 
J-b 
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Lemma 3.2 (A priori bound for the mass) For all a > 0, b > bo > 0, 1/2 < r < 1, any 

solution (c, u) of P T (a,b) satisfies 

f b i \ ( 2 maxR r \ . 

— J u ' x ' z ) — max I — — , m rj, vie [- a, a|. 

Proof. Integrating w.r.t. z the inequality < (maxR r — m(x)) w + (i? 2 + l) u zz — 

2Bu xz , we get 

—m"(x) — cm! (x) < ^maxr — — fc~ m(a;)^ m(x) + (-B 2 + l) (u z (x, b) — u z (x, —b)) 

-2B (u x (x, b) - u m (x, -b)) . 

Since u z {x, b) < 0, u z (x, —b) > and u x {x, b) = u x {x 1 —6) = 0, the mass satisfies the Fisher-KPP 
inequality —to" — cm! < (maxg r — to) to. Since m(— a) < top (see (19)) and m(a) = 0, the 
maximum principle concludes the proof of the lemma. □ 

The above nonlocal control now provides the following a priori bound for u. 

Lemma 3.3 (A priori bound for u) There exists M > such that, for all a > 0, b > bo > ; 

< t < 1, any solution (c, u) of P T {a 1 b) with < c < c* + 1 satisfies 

< u(x, z) < M, V(x, z) e [-a, a] x [-b, b] . 
Proof. For < r < 1/2, one keeps the local part and writes —£{u) — cu x < (maxRr — ^u) u; 



recalling (18), the maximum principle then gives a control of u by max ^ 2 m ^ x r j C 



^2(5/3 



Next, for 1/2 < r < 1, let us denote by (xm,zm) a point where u achieves its maximum M, 



and by B r the ball centered at (xm,zm) with radius r > 0. Note that < u < C 



^25/3 



on 



dQ. The function w := u — C 



^28/3 



therefore satisfies 



—£(w) — cw x — (maxi r)w < Co '■= (maxR r)C 
w < 



^28/3 



in Q 
on dQ. 



From the local maximum principle [25, Theorem 9.20] and it extension up to balls intersecting the 
boundary of the domain [25, Theorem 9.26], we infer that 



sup w < C\ 

B 1/2 nQ 



1 



B 1 nQ 



w + + C 2 \\Co\\ L »(B 1 nQ) 



where C\ = C\(B) and C2 = C2{B) are positive constants. Notice that C\ does not depend on 
c, which is a coefficient of the operator L(w) := £(w) + cw x + w, because c belongs to a bounded 
interval, namely [0, c* + 1]. Using successively u > and Lemma 3.2 we deduce that 



/ w+ - 1 


1 U + 


/ ° 


y2S/3 

00 


< 2 max ( 


1 BxC]Q J 


BiHQ 


J B±nQ 




00 \ 



' 2 maxR r 



m r + C 



p2(5/3 



I Si 



Recalling that M = maxu is achieved at the center of the ball B 1 / 2l we deduce form the upper 
estimates that 



M < C 



p2(5/3 



Si I 



max 



2 maxR r 



mr + C 



p2(5/3 



C2C0IS1 



This concludes the proof of the lemma. 



□ 



We now provide a control of the tails of the solutions as \z\ —> 00 by appropriate Gaussian 
functions (recall estimate (17)). 
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Lemma 3.4 (Gaussian control of the tails of u) There exists M > such that, for all a > 0. 
b > bo > ; < r < 1, any solution (c, u) of P T (a, b) with < c < c* + 1 satisfies 

< u(x, z) < MT 2 £ /3 {z), V(x, z) e [-a, a] x [-6, 6]. 

Proof. First observe that 

- £{u) - cu x - r{z)u < on Q, u < CT^ /3 on dQ. (22) 

Define <f>(x,z) = <p(z) :~ MT^ 3 (z), with M > to be specified later. Recall that solves 
(15) so that - c4> x - r(z)4> = (^-z 2 + A^ /3 ) 4>. Therefore, if (3 := max U, \j "^f ' ^j , we 

have 

- £(0) - c<fe - r(z)0 > on (-a, a) x (/3,b). (23) 

Let us now select 

M := max ( C, , M 2g/3 ) , (24) 
\ min^^j IV ) 

where M is as in the previous lemma. The choice (24) enforces 4>(x,z) > CT^ 3 (z) > u(x, z) on 
{=Fa} x [/3, b] U [—a, a] x {&}, and 0(x,z) > M > u(x,z) on [—a, a] x {(3}. Hence the comparison 
principle — note that r(z) < 0, when z > /3 — yields u < <f> on [—a, a] x [/3,6] and, by the choice 
(24), on [—a, a] x [0,6]. Similarly, u < (f> on [—a, a] x [—6,0]. The lemma is proved. □ 

3.4 A priori estimates for c 

We provide a priori bounds for the speed c of solutions to P T {a,b). We first show that, roughly 
speaking, too rapid waves solutions of P T (a,b) have too small value at (x, z) = (0,0). We recall 
that the speed c* was defined in (9). 

Lemma 3.5 (A priori upper bound for c) Let b > and e € (0, 1) be arbitrary. Then there 
exists do = ao(e, b) > such that, for all a > ao, all < r < 1, any solution (c, u) of P T (a, b) with 
c > c* satisfies u(0, 0) < e — and therefore cannot solve P(a,6, e). 

Proof. Let b > be given. Assume c > c* and let us show that u(0, 0) — > as a — > oo. 

The function u satisfies —£(u) — cu x — r{z)u < in Q. Therefore, changing variables, the 
function 

v(x, y) := u (^p=, yf& + \y) , (25) 

satisfies 

- v xx - v yy - c^B 2 + lv x - r (y/B* + ly) v<0, (26) 

■-By 



Qi ■= {(x,y) : \y\ 



< 



< a 



We shall now construct a positive solution. Since c > c* , one can select /j, < such that 
fi 2 + cy/B 2 + 1/i + ^ (B 2 + l) = and define <p(s) := e^ s which solves 

c* 2 



- c\jB 2 + lip' - — (B 2 + 1) <p = 0. 
Now, let us define 

r 5 / 3 / . \ 

n / / \ b ill aVB 2 + l — =B4= 

w(x, y) := Ka <p(x)r^ ( y/B 2 + ly) , n a := ^- e V vTO^ 

V / mmr >, m1" 
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with the eigenfunction appearing in Definition 1.2. Using direct computations and the defini- 
tion of c* in (9) we see that 



IVxx ^yy 



\/B 2 + iw x - r ( \/B 2 + 1 y ) w = in Q x . (27) 



We now compare the values of v and w on the boundary. If (x, y) € 9Qi is such that ~7ggj^ 7^ 



-a, then v(x,y) — < w(x,y). If (a;, y) e 9<2i is such that 7gg=j= = —a then 



= — /Hf tVlATl 



v(x,y) = T d b /6 UB 2 + ly) < K a e V vTO; min < w(x, y). 

\ / [— 6,6] 

As a result, we have v < w on 

Now a classical argument will imply »<ioon the whole of Q%. Indeed, since v is bounded and 
w > on Qi, we can define ao '■= max {a > : av < w in Qi } > 0. Then a$v < w and there is 
a point (xo,yo) such that aow^o, 2/o) = w ( x o,yo)- I n view of (26), (27) and the strong maximum 
principle the point (xq, yo) has to lie on dQi, which enforces ao > 1. Thus v < w in Q±, and 

u(0, 0) = v(0, 0) < w(0, 0) = K a ^0asa^oo, 

which concludes the proof of the lemma. □ 
Next, we show that standing waves (i.e. c = 0) have too large value at (x, z) = (0, 0). 

Lemma 3.6 (Standing waves: a priori lower bound for u(0, 0)) There is e* > such that 
if a, b are large enough, then, for all < r < 1, any standing solution (c = 0,u) of P T {a,b) 
satisfies m(0,0) > e* — and therefore cannot solve P(a,b,e) for any e € (0,e*). 

Proof. For R > 0, let us introduce (hr, Yr) as the solution of the principal eigenvalue problem 

-£(Y R )(x, z) - r(z)Y R {x, z) = fi R Y R (x, z) for all (x, z) e (-R, R) 2 
Y R = ond((-R,R) 2 ) (28) 
Y R (x,z)>0 for all (x,z) e (-R,R) 2 , T fl (0,0) = l. 

If R = oo, the above problem is equivalent to (8), and [i R —> A^, as R — > oo. Let us therefore fix 
R > large enough so that 

< m < < 0, (29) 

and 

MT^{z)dz<^. (30) 

R,R]° 4 

Next, let a > R + 1, b > R + 1, < r < 1 be given, and (c = 0, u) be a solution of 
P T (a,b). Thanks to the Harnack inequality, there exists C > (independent of a and 6), such 
that ||«||.L°°([-ii,iip) < Cu(0, 0), and then 



max 



k(z, z')u(x, z') dz' + 7 (1 -t)u) <C (2Rk+ + 7) u(0, 0) + ' 



(x,z)e[-R,R] a \ 7_ b ' I v ; 4 ' 

where we have used Lemma 3.4 and (30). As a result, u satisfies 

-S(u) - r{z)u > (-C (2Rk+ + 7) u(0, 0) + u in (-R, R) 2 , 

and u > on d ((-iZ, R) 2 ). Hence, if 
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u becomes a super-solution for (28) — thanks to (29). We conclude as in the proof of Lemma 3.5: 
defining ao := max {a > : aT r < u in [-R, R} 2 } > and using the strong maximum principle 
we see that aoTji = u on [— R, R} 2 , which contradicts u > 0. It follows that u(0, 0) > e* . The 
lemma is proved. □ 

3.5 On some related local problems 

First, we show the well-posedness of the local problem Po(a, b). 

Lemma 3.7 (Well-posedness for r = 0) There exists 7 > such that, if a, b are large enough 
and e £ (0, £*), there exists a unique (c, u > 0) £ R x W 2 '°°(Q), solution of Po(a,b), namely 



P (a,b) 



—E{u){x, z) — cu x (x, z) = (r(z) — 7 u(x, z)) u{x, z) in Q 
u(x,z) = l {x= _ a} (x)T S b /3 (z) on dQ, 

such that w(0, 0) — e. Moreover, by the above a priori estimates, < c < c* and < u < M . 



Proof. A standard argument proves that there is a unique positive solution to Po(a,b). For the 
convenience of the reader let us prove this fact. Since and a large enough positive constant are 
respectively a sub- and a super-solution of Po{a, b), the existence of a positive solution to Po{a, b) 
can be obtained using a classical monotone iterative scheme. Also, by the maximum principle, 
any positive solution of Po(a,b) is bounded. Now let u and v be two bounded positive solutions 
of Po(a, b). Thanks to the boundary condition and the Hopf lemma the following quantity is well 
defined: 

r* := inf {t > : V(a;, z) £ Q, u(x, z) < tv(x, z)} . 

Assume by contradiction that r* > 1. From the definition of r* , the boundary condition and the 
Hopf lemma, there exists (xq, zq) £ Q such that u(xq, zq) = t*v(xq, Zq). At this point, we get the 
contradiction 

< £(t*v - u)(x ,z ) + c(t*v - u) x (x ,z ) < JT*(1 - t*)v 2 (x ,z ) < 0. 

Thus t* < 1 and we have u < v. By interchanging the role of u and v, we get u = v. Hence there 
is a unique positive solution to Po(a, b). 

In order to apply a sliding method, we slightly modify the Dirichlet boundary conditions and 
consider, for small 77 > 0, 



PS(a,b) 



—£{v rt ){x, z) — cv r J,(x, z) = (r(z) — ^v n {x, z)) v n (x, z) in Q 
! ;"(, 1 z) = ^rf'"(z) ondQ, 



where T^ 3 ' v (z) :— T^ 3 (z) + •qTt^(z). Note that for any nonnegative rj, and a large enough 
positive constant are respectively a sub- and a super-solution of Pq (a, b). Thus the existence of 
a positive and bounded solution of Pq (a, b) can be obtained using a classical monotone iterative 
scheme. We shall now select 7 > so that (x, z) i-> r^ 3 ' v (z) become a super-solution for P^a, 6). 
Using (15), (16), X S b /3 > A^ 3 and r£ A '' > r^ /3 we see that 

- (b 2 + 1) A.rf^ - (r(*) - 7 i1 73 <") r s b ^> > ^ + S -z 2 + 1 r 5 b ^ r 5 b ^>, (31) 



which is clearly nonnegative for \z\ > z := y g . Now, as soon as b > z + 1, it follows from 

f 3 (z) > ^rf (o) = ^ 



the Harnack inequality that there is C > such that T^ 3 (z) > ir^ 3 (0) = i for all \z\ < z. 



5 /3 

Hence if we select 7 := —C\Jo > 0, the right-hand side member of (31) becomes nonnegative for 
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\z\ < z, and (x, z) i— > r^ 3 '''(z) is a super-solution for P$(a, b). Next, we show that, any solution 
of Pg (a, 6) satisfies 

<lt /3 '"(*), V(v)eQ. (32) 

Indeed, for w' 7 a non negative solution of Pq(o, 6), the following quantity is well defined 

a := sup{a > : V(x, z) € Q, av n (x,z) < rf /3,, '(z)} <= (0,1], 

since T S b /3 ' v > in Q. Let us assume by contradiction that cto < 1. In view of the boundary 
conditions for z; 17 , this implies that a point (xo,^o) where aoV n (xo, Zq) = r^ 3,, '(zo) cannot be on 
dQ. Hence w := r^ 3 ' v — aov n has a zero minimum at (xo, zq) € Q and 

> (-£ (u>) - cw x )(x ,z ) > ~fa (v v ) 2 (x a ,z a )(l - a ) > 0, 

which is absurd. Hence ao = 1, and < w'' < r^ 3 ' ?? in Q. Then (32) follows by applying the 
strong maximum principle. 

By the classical sliding method [11], v v is strictly decreasing in the x variable. For the con- 
venience of the reader, let us give a proof of this fact. For h > 0, define vV{x,z) := v v (x + h, z). 
Since v\ a = < v v thanks to (32), one can define 

h* := inf {h > : Vr g [/i, 2a], <' < w''} . 

Assume /i* > 0. Then there are sequences h n /* h* , {x n ,z n ) with (i„,z n ) > v r '(x n , z n ). 
After extraction and using that the infimum /i* is achieved, we have a point (xooi^oo) such that 
Vh* (xooj Zoo) = v^ixoo, Zoa), i.e. a point of zero maximum for u?, — v''. Because of the boundary 
conditions the point (xoo, Zqo) cannot lie on the upper or the lower boundary of (—a, a—h*)x (—b,b). 
In view of (32), it is neither allowed to lie on the left or right boundary of (—a, a— h*)x (— 6, 6). Since 
w r? and uJL are both solutions of— £(v) — cv x = [r(z) — jv) v in (— a, a— h*) X (— b, b), the maximum 
principle then yields i>?» = w 17 , i.e. v T, (x,z) = v v (x + h*,z). Hence v v (x, z) = v v (x + nh*,z) for 
all n > 0. Letting n — > oo yields v v = 0, which is a contradiction. It follows that h* = and v v is 
non increasing in the x variable. 

Now, we construct a solution to Po(a, b) as a limit, as 77 — > 0, of solutions to Pq(q, b). The 
interior elliptic estimates imply that, for all 1 < p < 00, the sequence (v v ) is bounded in W 2 ' P (Q). 
From Sobolev embedding theorem, one can extract a subsequence (u'') converging to some u, 
strongly in C 1 ^(Q) and weakly in W 2 ' P (Q). Moreover, u is a solution of Po(a, b). As a limit 
of decreasing functions, u is decreasing in the x variable. By differentiating the equation and 
applying the maximum principle, one then obtains the strict decreasing of u w.r.t. x. 

It is then standard that, if (ci,ui) and (02,1*2) are two solutions of P^{a,b) with c\ > C2, 
then Mi < U2- Indeed, ui is a super-solution of the equation for {c\,u{). Hence there exists a 
solution for this equation which is below ui. By uniqueness this solution is u\. Hence u\ < ui 
and, by the strong maximum principle, u\ < U2- As seen in Lemma 3.5, Lemma 3.6, the solution 
of Po(a,b) with speed c = 0, c > c* satisfy u(0, 0) > e* > e, u(0,0) < e respectively, if a, b are 
large enough. Then, there is a unique c, which belongs to (0,c*], such that the solution (c, u) of 
Po(a,b) is e- normalized. The lemma is proved. □ 

In order to apply a Leray-Schauder degree argument in the next subsection, we also need to 
consider the family < a < 1 of local problems 

!—£(u)(x, z) — cu x (x, z) = (r(z) — (1 — a)7Z — cryu(x, z)) u(x, z) in Q 
(33) 
u(x,z) = l {x= _ a} (x)T b ' {z) ondQ, 

where 1Z := max 2e u (r(z) + |z 2 ). 

Lemma 3.8 (On local problems P a (a,b)) (i) There exists M > such that, for all a > 0. 
b > 0, < a < 1, any solution (c, u > 0) of P a {a, b) satisfies < u < M. 



12 



(ii) Let b > and e G (0, 1) be arbitrary. Then there exists ao = Oq(e, b) > such that, for all 
a > ao, all < a < 1, any solution (c, u > 0) o/ P a (a, b) with c > c* satisfies u(0, 0) < e. 

(m) There exists £o > 0, ao > such that, for any a > ao, i/iere exists a speed — c = — c(a) < 
smc/i £/ia£ /or aH 6 > 1, a// < a < 1, any solution (c, u > 0) of P a (a, b) with c < — c satisfies 
u(0,0) > e - 

(iu) //a, & are ^arge enough and e G (0,£o)> i^en t/iere exists a unique (c, u > 0) € K x W^ 2, °°((5) 
solution of 

!—£(u)(x, z) — cu x (x, z) — (r(z) — lZ)u(x, z) in Q 
5/3 (34) 
u(x,z) = l {2;= _ a }(a;)r b / (z) on <9Q, 

wi£a it(0, 0) = e. Moreover, by the above a priori estimates, — c < c < c* and < u < M . 

Proof. Item (i) follows from —£(u) — cu x < a(lZ — ju)u, (18) and the maximum principle. 
Define 

»<*.») == - (7=^. ^») ■ * == {<*■»> = M < 78=. |7p=^| < ■} ■ 

Then v is a subsolution of (26) and, to prove (ii), we can reproduce the proof of Lemma 3.5. 
Let us prove (in). Observe that v solves 

— v xx — v yy — c\J B 2 + lv x = (r(\/ B 2 + 1 y) — (1 — cr)1Z — a^/vj v in Q\, (35) 

that v(x,y) = for (x,y) G dQi such that ^jpjjj^ 7^ —a, and that v(x,y) = T^ 3 (V B 2 + 1 y) for 
(x,y) G 9Qi such that ^gsjn = —a- It follows from (16) and the Harnack inequality that there 
exists C > such that 

L* /3 (z) > CLf /3 (0) = C, for all 6 > 1, all \z\ < 1. (36) 

Define, for a > 0, 

C- ^ x + a +(aVB J TT+ 7 # T7i 



max[_ M] I 2aVB 2 + 1 

inQ 2 :={(x,y): \y\ < 7 ^=, < a} = QiH {|y| < 7^}- Since t =2aVWI <0< « 

in Q2, we can define 

ao := min {a G 0, 2a \J B 2 + 1 : V(x, y) € Q2 , ip a (x, y) < v(x, . 

Assume by contradiction that ao > 0. Observe that ?p ao (x,y) < = v(x, y) for (x, y) G <9Q2 
such that ^jpjjj^ = a an d < ^ga +1 > tnat ^aoO^f) = < v(x,y) for (x,y) G 9Q2 such that 
|Z/| = 7^= and -a < < a, and that - by (36)- ^ ao (x,y) < T S b /3 (VW+Ty) = v(x,y) 

for (x, y) G 9Q2 such that Jg-Pj^ = — a - Therefore the only two points on 9Q2 where v — ijj ao may 
attain its zero minimum value are I . =k ± tM ± ) . But since ao > we see that, for 

e > small enough, ip ao - e < < f in a neighborhood of these two points of the boundary 9Q 2 - 
Hence, by the definition of ao, there must be a point (xo, yo) S Q2 where v — ip ao attains its zero 
minimum value, so that > —A x ,y ( v — V'qo) (^Oj Vo) — cV B 2 + ld x (v — ip ao ) ( x o, Ho)- Using (35), 
(15) and straightforward computations, we arrive at 

> v(x , y ) (-(1 - a)Tl — A? — a 1V (x , y )) - r° U B 2 + 1 y Q ) 



> vi x ,y )(-K-X°-,M)-±Cr°(VB^ly ) 



2a max 
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Since v(x ,y ) = ip ao (x ,y ) < ma ^ r u ^° (V B 2 + ly ) and since -Tl-\\ < -max[_ M j r-A? < 
we end up with 

^d^U^ 1 ") (- R -^ M -£)- 

Hence, for large negative speed, namely 

c < -c = -c(a) := 2a (-ft - A? - jM) < 0, 
we get a contradiction, so that «o = 0. It follows that 

u(o, o) = v(o, o) > Mo, o) = r0 U - irzjm , tt ) > — = : £ o, 



maxj.!,!] F° V2 2a(B 2 + l)/ 3max hlil] r0 

for a > ao, with ao > sufficiently large and independent on 6 > 1 and < c < 1. This concludes 
the proof of (Hi). 

To prove (iv), observe that, since —1Z < A^ 3 , r^ 3,r '(z) :— T^ 3 (z) + r]Tti 3 (z) is a supersolution 
for (35). Hence we can reproduce the proof of Lemma 3.7. Notice in particular that if (ci, u\) and 
(02,1*2) are two solutions of Po(a,b) with c\ > C2, then ui < u^. □ 

3.6 Construction of a solution in the box 

Equipped with a priori estimates of Subsections 3.4 and 3.3, we are now in the position to construct 
a solution to P(a, b, e), with e G (0, min(e*, £o)). We shall use a Leray-Schauder topological degree 
argument (see e.g. [10] or [9] for related arguments). 

Proposition 3.9 (The solution in a box) Let e £ (0, min(e*, £o)) be arbitrary. There exist 
K > and bo > such that for any b > bo the following holds. There exists a® = ao(b,e) such 
that, for all a > ao, the problem P(a, &, e) has a solution (c, u) such that 

\\u\\ C 2 {Q) < K, 0<c<c*. (37) 

Proof. For a given nonnegative function v defined on Q and satisfying the Dirichlet boundary 
conditions as requested in P(a, b, e), consider the family < r < 1 of linear problems 



(PI) 



—£(U)(x, z) — cl/ x (x, z) 

= [ r(z) — t I fc(z, z')v(x, z') dz' — 7(1 — t)v(x, z) J v(x, z) in Q ^gg-j 



-fa 

5/3 



lU(x,z) = l {x= - a} (x)r° b '*(z) ondQ. 

Let us define IC T the solution operator of the above system. More precisely JC T is the mapping of 
the Banach space X := M x C 1,a (Q) — equipped with the norm |j(c, v)\\x '■= max(|c|, HvH^i.c) — 
onto itself defined by 

fCr ■ (c, v) ^ (e - v(0, 0) + c, Ul := the solution of (P C T )) . 

Constructing a solution of P(a, b, e) is equivalent to showing that the kernel of Id — /Ci is nontrivial. 
The operator IC T is compact and depends continuously on the parameter < r < 1. Thus the 
Leray-Schauder topological argument can be applied. Define the open set 

S := {(c,v) : < c < c* + 1, v > 0, \\v\\c^ < M + 1} C X, 

where M > is as in Lemma 3.3. It follows from the a priori estimates Lemma 3.3, Lemma 3.5 
and Lemma 3.6, that there exists ao = ao(b,e) > such that, for any a > ao, < r < 1 the 
operator Id — JC T cannot vanish on the boundary dS. By the homotopy invariance of the degree we 
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thus have deg(Id — Ki,S,0) = deg(Id — /Co, S, 0). Additionally, thanks to Lemma 3.7, any element 
of the kernel of Id — /Co belongs to S, so that 

deg(Id - /Co, S, 0) = deg(Id - /C , S, 0), 

where 

S := {(c,v) : -c < c < c* + 1, v > 0, \\v\\ C i,a < M + 1} C A, 

with c ~ c(a) > as in Lemma 3.8 (Hi). 

Let us now consider the family < a < 1 of local and linear problems associated with (33), 
namely 



(Pc) 



—£(U)(x, z) — cU x (x 1 z) = (r(z) — (1 — a)lZ — ajv(x, z)) v(x, z) in Q 
U(x, z) = l {x= _ a} (x)T S b /3 (z) on dQ, 

and let K. a be the associated solution operator, that is 

t a : (c, v) ^ (e - v(0, 0) + c, W c := the solution of (P a c ) 



The operator K,„ is compact and depends continuously on the parameter < a < 1. The analysis 
of the local problems P a (a, b) in Lemma 3.8 shows that Id — IC a cannot vanish on the boundary 
of S. Since /Co = ICi we have 

deg fid - /Co, S, O] = deg fid - fc Q , S, 



To complete the proof, let us compute deg (Id — /Co, S, 0J by using two additional homotopies. 
First, consider, for < r < 1, 

g T : (c, v)^(e-(l- t)v(0, 0) - tU°(0, 0) + c, U° := the solution of (P c )) . 

If Q T (c, v) = (c,v) for some (c, v) £ dS, then (c, v) £ dS solves the local problem P°(a,b) and 
is such that v(0, 0) = e. By the a priori estimates of Lemma 3.8, this cannot be. Therefore 

Id — Q T does not vanish on the boundary dS. Since X^o = Go we have deg fid — K.q, S, 0^ = 
deg f^Id — Qi, S, 0^. Next, we know from Lemma 3.8 (iv) that there is a unique (co, U Ca ) G S which 
solves the local problem P°(a, b) and is such that U Co (0, 0) = e. Then, consider, for < t < 1, 

n T : (c, v) h> (e - ?7 r (0, 0) + c, r?7 c + (1 - r)?7 Co ) . 

If 7i T (c,v) = (c,v) for some (c, u) € 95*, then the uniqueness in Lemma 3.8 (iv) enforces c = Co, 
£f? = U CQ = v so that (c, w) e 9S 1 solves the local problem P°(a,b) and is such that v(0, 0) = e, 
which cannot be. Therefore Id — H T does not vanish on the boundary dS. Since Hi = Gi we have 

deg (id - Q x , S, o) = deg (id - 5, o), where 

Id - T^o : (c, «) ^ (f/ c °(0, 0) - e, w - . 

As seen in the proof of Lemma 3.8 (iv), C/°(0, 0) is strictly decreasing in c so the degree of the 
first component of the above operator is —1. Clearly the degree of the second one is 1. Hence 
deg(Id — Hq, S,0) = —1 so that deg(Id — JCi,S, 0) = —1 and there is a (c, u) € S solution of 
P(a, b, e). This concludes the proof of the proposition. □ 
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4 The front with minimal speed c* 



Equipped with the solution (c, u) of P(a, b, e), with e € (0, min(e*, £o)); of Proposition 3.9, we now 
let a —> +oo. Note that we have the bounds (37) on c and u, and also the Gaussian control of the 
tails in Lemma 3.4. This enables to construct — passing to a subsequence a„ — > +oo — a speed 
< Cb < c* and a function Ub £ C t 2 (1 x [—6, b]) with the same bounds as those of u. Similarly, 
we can then consider b — > +oo, to construct, via a subsequence b n — > +oo, a speed < c < c* and 
a function u g C 2 (R 2 ), such that < u < K, and 

— £(u)(x, z) — cu x (x, z) = I r(z) — / k(z, z')u(x, z') dz' ) u(x, z) in R 2 (39) 



u(0,0)=£ (40) 
< u{x,z) < Mrf(z), V(x,z)eM 2 . (41) 

4.1 The constructed wave has the minimal speed c* 

Here, we show that, by reducing the normalization (40) if necessary, the above constructed solution 
has speed c = c* . 

Lemma 4.1 (A priori estimate for the infimum) There exists e > such that any solution 
(c,u) of (39) ; (41) with c > and inf^gg u{x 1 0) > actually satisfies inf^R u(x, 0) > e. 

Proof. We choose R > large enough, such that 

A^o < A° fl < % < 0, fc+ / Afr^/ 3 (z) cfo < (42) 

Thanks to the Harnack inequality, there exists C > such that 

< Cu(x, z) for all x e M, |z| < R, \z'\ < R, (43) 

— A 

which, combined with (41) and the second part of (42), implies J k(z, z')u(x, z') dz' < 4 " + 
2k+CRu(x,z) in the strip R x (-R, R). Hence, -S(u) - cu x > (r{z) - (— ^ + 2k+CRu)) u in 



x 



(— R, R). Changing variables, the function v(x, y) := u ^-^=J=K=, \JB 2 + 1 y^j then satisfies 



-A° 



v / S 2 + lu a: > rfV^ + ly)- — -^ + 2fc + C*i?w u, (44) 



Now let 77 > be arbitrarily given. Define, for a > 0, 



^ a (a,y) := «(1 - V^r (\AB 2 + ly) 



Observe that the Harnack inequality (43) implies ip a <v for a = C 1 ||r5j|| inf^gg u(x, 0), and 
that (41) implies ipa(0,0) = M > u(0,0) for a = M. We can therefore define 

«o := max {a > : V(x, y) S 5, ip a (x, y) < v(x, y)} & (0, M]. 

Hence v — i\i ao attains a zero minimum at a point {xq, yo) — depending on n — which must lie in 



V-B 2 + l ' Vb^+T 



so that 



> -&x,y (v - ^ao) (zo,Z/o) - c\/ S 2 + ld x (v - ip ao ) (x ,y Q ). 
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Using (44) and straightforward computations, we arrive at 
> 



\0 

_ ih+r 



- 2k + CRv(x ,y ) - \°r) v(x , yo) - 2a (r? + c^J B 2 + lryxo) T° (V^ 2 + lj/ ) . 
Using the first part of (42), «o < M an d |xq| < this yields 

0> (^-2k + CRv(x ,y )^v(x ,y )-2M\\T R \\ oo (r] + cy/B* + 10?) • 
It follows from the Harnack inequality (43) that w(xo,yo) > ■h inf^gH m(x, 0) > 0, so that 



nz ,yo) > ou+r ,r> — ; ^ — r~?v; • 

2k+CR y 4 mf^gR m(x, 0) / 

Since ?] > can be chosen arbitrarily small, we have v(xo,yo) > sk+CR anc ^ then ao > e 



8i+CR r» 



> 0. Hence, u(x,y) > e(l - r)x 2 )T° R [\l B 2 + ly) for all (x,y) £ S. Since rj > can 



be chosen arbitrarily small, we have v(x, y) > §r^ (V _B 2 + ly) , and in particular inf^gg u(x, 0) = 
inf-ugR v(x, 0) > |. This proves the lemma. □ 

As a result, the constructed solution of (39), (40), (41) satisfies inf l6 i u(x, 0) = 0. Without 
loss of generality, we may assume lim inf K _>oo u(x, 0) = 0. The following proposition then enforces 
c = c* for the constructed wave. It is also of independent interest since it proves the non existence 
of waves for < c < c* as stated in Theorem 1.4 (ii). 

Proposition 4.2 (c = c* for the constructed wave) Any solution (c,u) of (39), (41) withe > 
and lim inf ^-j.oo u(x, 0) = actually satisfies c > c* . 

Proof. Assume by contradiction that < c < c* . Choose c < c < c* . Since X R — > A^ = 
— {B 2 + 1)%- as R — > co, we can choose R > such that 



*i_ < _£_i^ i k+f STT< &±±i^jt. (45) 

B 2 + 1 - 4 2 4 ' ' 



-R,R]<= 



Let us define the open rectangle 



Thanks to the Harnack inequality, there exists C > such that for any solution (c, u) of (39), (41) 
with < c < c*, and for all Xi € M 

max u ( xi + ^=£L, z' ) < Cn(xi, 0). (46) 



(avy,z')efix[-ii,ii] V V-B 2 + 1 

Following the change of variables of Lemma 3.5, we see that v(x, y) := u ( , V B' 2 + 1 



satisfies 



cV-B 2 + lw x 



(V^ 2 + ly) - /" fc(Vs 2 + ly,z') 



a: - By f . 



> 



B 2 + 1 c* 2 - 5 2 



VB 2 + ly) - ^ — - + k+2CRu(0,0) 
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where we have used the second inequality in (45) and (46) with x\ = 0. Next, define the function 

,/ x n / r-z \ c^Tl* ( J(B 2 + 1) (d 2 - c 2 ) tt\ 

i>(x, y) := T% (y& + ly) e — sin I ^ ^ '-x + - j . 

We have ip = on dil and, using the first inequality in (45), 

, N / ^2 \ £? 2 + 1 c* 2 — c 2 

-VW - V«w - cV-B 2 + - r (yW+ly) V = i(B 2 + 1)- + A^J V < 2 4 V>, 

in 0. 

Since < v < K, we can define ao := max{a > : aip < v in f2} > and there is a point 
(^Oi^o) € where w := aoV' — v attains a zero maximum. In view of the above inequalities, we 
have at point (xo,yo), 



< -W xx - w yy - c\J B 2 + lw x - r (y/ B 2 + 1 y j w 

B 2 + 1 tf 



fc+2Ci?u(0,0) 



v(x ,yo), (47) 



which in turn implies w(0,0) > (B 2 + 1) ^k+CR ' ^ n v ^ ew 01 (46), the argument is invariant under 
translation w.r.t. x variable, so that 

inf^0)> (B 2 + l)|^>0, 
which contradicts lim inf x ^oo u(x, 0) = 0. The proposition is proved. □ 

4.2 Behaviors as x — > ±oo 

The following proposition will show that the constructed wave satisfies the lower bound in (11) 
and (12) in Theorem 1.4. 

Proposition 4.3 (Behaviors at infinity) Let (c, u) be a solution of (39), (40), (41) with c > 0. 
Then the following holds. 

(i) There exist R > and k > such that u(x,z) > kT^z) for all (x,z) G (— oo,0] x [—R,R]. 

(ii) If e > is as in Lemma 4-1, then both I u{x, z) dz and max zg RM(x, z) — > 0, as x — > oo. 

Jn 

Proof. Let us prove (i). We start as in the proof of Lemma 4.1: choose R > large enough so 
that (42) holds, choose C > such that both (43) and 

1 e 

min v(x, y) > -u(0, 0) = - (48) 

-i< x < , Ivl^^. C C 

hold, and observe that v(x,y) := u ^ J-g^x i V-B 2 + 1 2/) satisfies (44) in the strip S := R x 
V^+T ' v^+i) ' Now for r > > °' we define 

1 \ „n / r^T. T \ ( £ — A„ 



^(^i 2/) : = a Q + ^ r R (V-B 2 + 12/) j a:=min^ 



cmwjsk+cRmw^ 



The definition of a then enforces ipi < v on S-, where S- := (— oo,0) x ^— , j ■ We 

can therefore define 

770 := min {n > : V(x,y) £ S- , ip v {x,y) < v(x,y)} G [0, 1]. 
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Let us assume by contradiction that 770 > 0. Function v — tp ri0 then attains a zero minimum at a 
point (xq, yo); the definition of a and the Harnack inequality (48) prevents xq = so that (xq, yo) 

has to lie in (-^.OJ x We therefore have > -A x , y (w - ipr, )(x ,y ) - 

c\/ -B 2 + ld x (v — ij) Vo ) (%o i Uo ) • Using (44), 8 x i(j Vo > and the first part of (42) we arrive at 

> ^ _ 2k + CRv(x ,y )\ v{x ,y ), 

which in turn implies 8k +^ R < ^(^Ojyo) — ipri (xo,yo) < § ll^fill ' w hich contradicts the 
definition of a. Hence 770 = and v(x,y) > (VB 2 + 1 y) for all (.T,y) 6 (— 00, 0] x 

— — r=§q=p . This concludes the proof of (i). 

Thanks to the Harnack inequality and the control of the tails u(x, z) < Aff^ 3 (z), in order 
to prove (ii) it is enough to prove u(x, 0) — >, as x — > 00. Assume by contradiction that there 
exists v > and x n — > +00 such that u(x n ,0) > u, for all n. Then, the proof of (i) shows 

that u(x, 0) > I min f y 1 5k+CR\\r° |[ ) ^ or a ^ x ^ (^iin)i and then for all x G R. 
Hence inf^gR u(x, 0) > so that Lemma 4.1 implies inf x ^_ u(x, 0) > e, which contradicts the 
normalization (40). This proves (ii). □ 

5 Faster fronts (c > c*) 



In this section we fix c > c* = 2y -gsqrj- an d construct a nonnegative function u € C (R ) solution 
of 

— £(u)(x, z) — cu x (x, z) = I r(z) — I k(z,z')u(x,z')dz' ) u(x,z) in K 2 . (49) 



Using the change of variables v(x, y) := u ( ^ B ?^ 1 v 7 B 2 + ly^j , we need to construct a nonnegative 
f = v(x,y) solution of 

Cv(x,y) := -v xx (x,y) - v yy (x,y) - c\J B 2 + lv x (x,y) - r \ J B 2 + ly) v(x,y) = 

- v(x, y) J^k(^B 2 + l y, z') v (x - By + -JL= Z ' , ^p=) dz ' in r2 - ( 50 ) 

Note also that solving the problem in the box 

— £(u)(x, z)— cu x (x, z) = \r(z)— I k(z,z')u(x,z')dz' u(x,z) in Q = (— a,a)x(— b, 6), (51) 



is equivalent to solving 



Cv(x, y) = -v(x, y) J k (V ' B 2 + 1 y, z'^j v (x - By + -j=^== z' , ^ = = ^ 



with Qi = 



0i = {(*•») 5 lvl<7fe,|^&|<«} 



dz' in Qi, 

(52) 



We first adapt the strategy of [9]: we construct a solution in the box by using sub and super- 
solutions and the Schauder fixed point theorem. This will allow to let a —> 00 but we shall need a 
extra argument to let b — > 00. 

Construction of sub and supersolutions. We use again the supersolution of Lemma 3.5: since 
c > c*, one can select \i < the largest root of + cV B 2 + lfi + ^j- (B 2 + l) =0. Then the 
function 

w(x,y) := e" 1 ^ ( VB 2 + ly) , 
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with the eigenfunction appearing in Definition 1.2, satisfies Cw(x,y) = in R 2 . 

Next, we aim at constructing a kind of subsolution. Precisely we look after a function h such 
that (note that the supersolution w appears in the integral term) 

Ch{x,y)<-h{x,y) J k(yW+ly,^w(x-By+ -^= z', ^rp[ ) dz ' in {h > 0} " 

(53) 

„ a z < 

Since fc < fc+ and / R e Vb 2 +i r^(z') dz' < oo, there is C > such that 

jf fc (>/F7l y, „ (x - fly + ^^z', ^^=) dz' < Ce^\ (54) 

Let us choose e > small enough so that — p := (p — e) 2 + c\/B 2 + l(p — e) + ^- (B 2 + l) < 
and p + e < 0. For a constant A > 1 to be selected later, let us define 

h(x,y) := (le"* - e^" e >) L^ (V^+Ty) = jw(x,y) - e^ x Tl (^W+ly) , 

with is nonnegative if and only if x > e _1 In. A. Thanks to the estimate (54), we have 

Ch(x, y) + h(x, y) J k [y/ B 2 + ly, z'^j w (x - By + -jJ^= z' , = dz' 

< (-p+Ce^ x ~ By) ^ h{x,y) 

< 0, (55) 

in {(x, y) e Qi : x > ~ In (£) + -^= b| which contains {/i > 0} = {(x, y) e Qi : £ > In A} 
provided that A is chosen sufficiently large. Notice that A does not depend on a but does depend 
on b. We will thus need an extra argument below. 

Construction of a propagating wave in the strip Ex [—6, b]. Let b > be arbitrary. Consider 
the problem (52) in the box Q\ supplemented with the boundary conditions 

v(x, y) = h<j(x, y) := max (0, h(x, y)) for all (x,y) <G dQ\. (56) 

Define the convex set of functions 

R a ,b :={wGC(Qi): h < v < w} , 

and the compact application & a ,b that maps a given v* £ C(Q±) to the solution v of 

£y(x, y) = -v(x, y) J ^ fc (V # 2 + ly, z') fx - By + -jJL= z ', _^==^ dz' i n Q l5 

supplemented with (56). Since £w < = Cw in Qi, i> = ho < w on dQ±, the maximum principle 
implies v < w in Qi. Also, (55) implies Ch$ + ho J kw < = £u + v J kv* < £v + v J kw 
in Qi H {ft, > 0}, and w > implies ^ < o on 9(Qi H {ft > 0}). Hence we have ho < v in 
Qi n {h > 0} and thus in Q\. Hence, <& a ,b maps R a .b into itself. By the Schauder fixed point 
theorem, <& a ,b has a fixed point v a ,6 € -R a ,& which solves the problem in the box (52) and satisfies 
the boundary conditions (56). Hence we are equipped with u a ^ solution of (51) with 



and 



u a , b (x, z) = h [^Jb 2 + 1 x + \~/W+t) ° n<9( ^' 



B z \ , n {-/W+1X+ 



(57) 



h i^VB^x + ^==z,^==j <u a ^ b (x,z)<e y Jrl(z) in Q. 

(58) 
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We claim (see proof below) that there exists M > such that, for all a > 0, b > 0, 

u a , b (x,z)<Mr 2 ^ 3 (z)<M T 2 ^ 3 =:M, in Q. (59) 

oo 

Now, for a given b > 0, choose A = A b as in the construction of h above. The family (u a _b) a 
is uniformly bounded, and is then uniformly bounded in C 2,a (Q). This allows to let a — > oo, 
possibly along a subsequence. In this limit, we have u a ^ — > Ub, which is a solution of (49) in the 
strip Sb ■= K x (—b,b), and satisfies (58), (59) in Sb', in particular (58) yields 

u b (x b ,0) > h Qln (j^-^Abjj =: e b > 0, x b := -^p = = ^ In (^"^f> 

Since the problem in the strip is invariant w.r.t. translation in the x variable, we may assume 
Ub(0, 0) > Eb > 0. Since e& — > as b — > oo, before letting b — > oo we need an additional argument 
to get a uniform w.r.t. b lower bound for w;,(0,0). 

Let us now prove (59). For < x < a, it follows from (58) that the mass satisfies m(x) := 

I-b u a,b{x, z) dz < C :— J R e v /b2+1 T^ (z)dz. For —a < x < 0, observe that w a .b = /io = on 
d((—a, 0) x (—6, &)) so that we can follow Lemma 3.2 to obtain that the mass satisfies a Fisher-KPP 
inequality. Since m(— a) = and m(0) < C, the maximum principle yields m(x) < max ( 2 ™^ - , C) 
for —a < x < and thus for —a < x < a. This uniform bound for the mass enables to argue 
exactly as in subsection 3.3 — recall that c > c* has been fixed — to get (59). 

Uniform lower bound for Uf,(0, 0). We choose R > large enough so that (42) holds, and 

C > such that (48) holds. Then, for b > R+ 1, v b {x, y) := u b ( ' v ' B2 + 1 y ) satisfies ( 44 ) 



in the strip S — M. x I — , , *l . We are therefore in the position to reproduce the proof 

\ V" +1 vB +1 / 

of Proposition 4.3 (z). Hence, there is re^ > 0, depending on Ub(0, 0), such that Vb(x, y) > KbT^y) 
for all (x, y) £ (— oo, 0] x [— i?, i?]. Since 6 > R + 1, we can apply the Harnack inequality to show 
that there exists Eb > such that 

Vb(x,y) > e b , V(x,y) S (-oo,0] x [-R,R], (60) 

Let us now introduce a smooth function <j> : R — > [0, 1], such that <f> = on (— oo, —4] U [—1, oo), 
and = 1 on [—3, —2]. For 77 > 0, a > 0, let us define 

i>a,r,(x,y) ■■= a(j)(r)x)T° R (V B 2 + lyj for (x,y) G 5. 

In view of (59) we have Vb < M so that, for any 77 > 0, we can define 

q = a {v,b) := max {a > : V(x, y) E S, ip a , v (x,y) < v b (x,y)} € [0,a max ], 

with a max := M/Ur^U^. Hence Vb — "0a o .»j attains a zero minimum at a point (xo,yo) ~ 
depending on r\ and — which must lie in (— ^, — ~) x ^— -^=|==, -^=|== J ; thanks to the definition 
of ^a.r; and 4>- Then, 



> -Ax,y (Vb - 1pa ,v) { x O,Vo) - C\/ B 2 + ld x (v b - 1p ao ,ri) (xo,Uo)- 

Using (44), the first part of (42) and straightforward computations, we arrive at 
> f— ^ - 2k+CRv b (x a , y )\ v b (x a ,y a ) + 



Q!0 



Then (60) yields 



> I — 2k^CRv b {x ,y ) — I v b {x ,y ) 
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Since 77 > can be arbitrarily small, we discover that Vb(xo,yo) > Sk +Q R , which in turn implies 
ao > e := 8k + c ^^u y — • Since ip ao , v (~§>o) = a ° there is a point Xb where Vb(xb,0) > e, which 
in turn provides a point x' b where Ub{x' b , 0) > e. Since the problem in the strip is invariant w.r.t. 
translation in the x variable, we can assume without loss of generality that x' b = 0. Thus we have 
the desired lower uniform bound Ub(0, 0) > e. 

Conclusion. The family {ub)b is uniformly bounded, and is then uniformly bounded in C 2,a (R 2 ), 
and we may pass to the limit b — > 00, possibly along a subsequence. In this limit, we have Ub — > u, 
which is a solution of (49) in K 2 , and satisfies (58), (59) in R 2 and u(0, 0) > e > 0. Hence, we 
have constructed u which satisfies (10), (13) — which in turn implies (12) — and the upper bound 
in (11). Last, the lower bound in (11) follows from Proposition 4.3 (i). This concludes the proof 
of Theorem 1.4 in the case c > c* . □ 
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